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J> Abstract 

l« The stable and unstable manifolds of a saddle fixed point (SEP) of the 

-^ Bonhoeffer-van der Pol oscillator are numerically studied. A correspon- 

rf\ dence between the existence of homoclinic tangencies (which are related 

to the creation or destruction of Smale horseshoes) and the chaos observed 
in the bifurcation diagram is described. It is observed that in the non- 
chaotic zones of the bifurcation diagram, there may or may not be Smale 
horseshoes, but there are no homoclinic tangencies. 

Keywords: BonhocfFcr-van der Pol oscillator; Smale horseshoes; chaos; bifur- 
cation; invariant manifolds. 
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1 Introduction 



The Bonhoeffer van der Pol oscillator (BvP) is the non-autonomous planar 
system 

x' = X- — -y + I{t) I ^ ^^^ 

v' = c.{x + a — hy) \ 

being a, 6, c real parameters, and / (t) an external forcement. We shall consider 
only a periodic forcement I{t) — A cos (27ri) and the specific values for the 




Figure 1: Bifurcation diagram for the first coordinate of the periodic points of 
the Poincare map, and 0.5 < A < 1.8. The blue Hne represents the position of 
the Saddle Fixed Point (SFP). 

parameters a = 0.7, b = 0.8, c = 0.1. These values were considered in [T] 
because of their physical and biological importance (see [2] ) . 

In previous works [3] , the existence of "horseshoe chaos" in BvP was studied 
analitically by means of the Melnikov method applied to an equivalent system, 
the Dufhng-van der Pol oscillator (DvP). It was concluded that the BvP system 
should have Smale horseshoes. Nevertheless, the method used there can be 
applied only for b > 1 and c < 1/6; which is not our case. In this work we will 
show the relation between the chaos transitions in the BvP system (fTl) and the 
creation or destruction of Smale horseshoes. Such horseshoes will be identified 
by the existence of homoclinic tangencies between the invariant manifolds of a 
saddle fixed point of the Poincare map. With this aim we perform a numerical 
descriptive analysis of the stable and unstable manifolds. 



2 The bifurcation diagram 

Obtaining, for different values of the parameter A, the periodic fixed points of 
the Poincare map (which is defined by the flow of the system on t = 1, see [4]), 
a typical bifurcation diagram is found, with chaotic and non chaotic zones (see 
Fig. [l]). Such diagram has been deeply studied in [5] and [T]. 

For ^ A ^ 0.61 and for 1.72 ^ A, there is a unique attracting fixed point. 
The first bifurcation takes place in A « 0.61, where two attracting 2-periodic 
points appear, and a saddle fixed point (SFP) between them. Our aim is to 
study the invariant manifolds of such SFP. 
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Figure 2: Some details of the bifurcation diagram of the x-coordinate. 

For 0.61 ^ A ^ 0.735 the bifurcation diagram is very simple, from a dynam- 
ical point of view. The attracting periodic points keep doubling its periods and 
forming a typical non-chaotic bifurcation diagram. Their attraction basins are 
perfectly separated and no fractal structures are formed. 

For 0.735 ^ A ^ 1.2835 there is a sequence of chaotic and non-chaotic zones 
(see Fig. [2]). The largest non-chaotic zone yields for 0.782 ^ A ^ 1.092, and 
there is a 4-periodic attracting point. For example, the four periodic points 
in the case A = 0.85 are xi w (-1.6444,0.4723), X2 ~ (-1.1557,-0.2346), 
X3 « (-0.8362,-0.4906), X4 w (1.8481,0.4416), and the image of any of them 
is the next one, i.e. f{xi mod 4) = Xt+i ,„od 4- 

Another event that must be mentioned is the sudden expansion of the attrac- 
tor that takes place at A ~ 0.748; once the horseshoe chaos has begun. Shuch 
expansion has been deeply studied in [T] using dynamical structure functions. 

For 1.2835 ^ A ^ 1.72 there is no horseshoe chaos. In this zone there are 
attracting periodic points. This situation is held until A « 1.72, where the SFP 
vanishes and a unique attracting fixed point appears. 

3 Invariant manifolds 



As we have mentioned above, our aim is to relate the chaos transitions in the 
bifurcation diagram to the creation and destruction of Smale horseshoes, which 
we shall identify with the existence of homoclinic tangencies between the invari- 
ant manifolds. In this section we are going to set the basic definitions and the 
notation related to the invariant manifods; then we shall describe the structure 
of the invariant manifolds, depending on the value of the parameter A. 




Figure 3: Attraction basins of the four attracting periodic points for A — 0.85. 
For such A, the SFP hes approximately in (—0.9384, —0.4146), which is in the 
boundary of the four attraction basins. 



Let xq be a SFP of a discrete dynamical system with a two-dimensional 
state space, given by a continuous map / : M^ — > M^. The stable manifold of 
xq, W (xo), is defined as 



W 



M :- { 



X e 



lim r{x)=xA 

n— >-oo J 



On the other hand, the unstable manifold of xq, W^ (xq)^ is defined as 



W^ (xo) := { 



X e 



lim r^{x) = xo} 

n— >-oo J 



In this case, these invariant manifolds are one-dimensional. Moreover, it is 
obvious that if they intersect at one point different from the SFP, then they 
must intersect at an infinite set of points and Smale horseshoes arc formed (see 

a)- 

A simple consequence of these definitions is that, in presence of several at- 
tracting fixed points, W (xq) is included in the boundary of the attraction 
basins of these points. Thus, if there are horseshoes, the attraction basins 
present a fractal structure (see Fig. [3|. 

Each invariant manifold has two branches that "arrive" at the SFP (in the 
case of the stable manifold) or "leave" the SFP (for the unstable manifold). 
More precisely, let x € W (xq), then d^ {/"■ {x) , xq) — >■ as n — > oo, being d* 
the arclength distance defined on W^ (xq). On the other hand, if a; G W" (xq), 
then d"(/"'(a;),a;o) — ?> oo, being d" the arclength distance defined on W^ (xq) 
(see El). 

In our case, the vector field / is given by the Poincare map defined by the 
flow of the system (fTj) on t = 1. This map has a unique SFP which appears 
with the first bifurcation at A w 0.61 and vanishes with the last bifurcation, 





Figure 4: Stable manifold for A — 0.70. Left: branches of the manifold plotted 
with the inverse method, before they get out of range. Right: stable manifold 
plotted with the method of "convergent points". The SFP is marked with a 
black circle. 

approximately at A « 1.72 (see Fig. fl]). Therefore we are going to focus on the 
values 0.61 ^ A ^ 1.72, which are the ones for which the invariant manifolds 
exist. 

3.1 The stable manifold 

The simplest method for plotting the stable manifold is the well known "inverse 
method" , which consists in plotting the unstable manifold of the inverse of the 
Poincare map. 

With this method, taking into account the machine precision we are using, 
we obtain that the branches of the stable manifold reach a numerical infinity at 
a finite time. This "out of range" takes place only for the x coordinate, due to 
the cubic term of ([I]). Nevertheless it is important to remark that, in theory, no 
point on the stable manifold can go to infinity at a finite time. 

Because of this behaviour, it is not possible to draw numerically the whole 
manifold using the inverse method. In fact, it only works until the branches are 
out of range (Fig. Hleft). 

Recently, new algorithms for plotting stable manifolds without computing 
the inverse have been developed (see [71 |B] ) . However this methods cannot 
be used here because they require that the manifolds are sufficiently bounded. 

To avoid these problems, in order to plot the stable manifold, we use another 
complementary method which consists in finding the points that, after a given 
number of iterations, lie near the SFP (Fig. H right). Nevertheless this method 
is computationally more expensive. For more details refer to Section [5J 



3.2 The unstable manifold 

The unstable manifold remains in a bounded region close to the SFP. However, 
the branches of the manifold present some folds that, at first sight, look like 
vertices, difficulting its study. 
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Figure 5: Unstable manifold for A = 0.70. The SFP is marked with a black 
circle. 



For example, for A = 0.70, the SFP is at (-0.9635,-0.4841), approxi- 
mately. One branch arrives to the zone near (—1.0704, —0.3385) where it folds 
(fold zone A). Other fold zones that have to be mentioned are near the points 
(-1.0545, -0.3848) (fold zone B) and (-1.0701, -0.341) (fold zone C). Since our 
Poincare map is orientation reversing, the other branch of the manifold has an 
analogous structure, presenting the fold zones A', B', and C, as the images of 
the fold zones A, B, and C. The branches are "trapped" between the fold zones 
B-C and B'-C respectively (see Figures [sj [6] and [t]) . 

If the parameter A is increased, the fold zones B, C, B', C approach the 
SFP (Figures |8J [o]) . When these fold zones arrive near the SFP, the unstable 
and stable manifolds will intersect, causing the first set of Smale horseshoes. 
For example, ii A = 0.735, these four fold zones are very close to the SFP (Fig. 
10). For A = 0.74 the first set of horseshoes has been created (Fig. 11). 



As we have mentioned in Section[2J for A « 0.748 there is a sudden expansion 
of the size of the attractor, and consequently of the size of the unstable manifold. 
For example, this expansion is well described in Fig. [Tljfor A = 0.75. In this 
case the main fold zones are near the points (1.7458, 0.2452), (-1.6781, 0.4752), 
(-1.2212,-0.2360), and (-0.8721,-0.5232), which are the zones of the bifur- 
cation diagram with greater density of points. 

In the non-chaotic big zone 0.782 ^ A ^ 1.092 there are four attracting 
points of period 4. These points are near the main fold zones as we can see for 
A = 0.85 in Fig. [TT} 



4 Smale horseshoes 

For 0.61 ^ A ^ 0.735 there are no intersections between the stable and unstable 



manifolds. So there is no horseshoe chaos in this region (Fig. 12). 

For A = 0.735 it is observed that several homoclinic tangencies are about to 
be formed. This will cause the creation of the first set of horseshoes and the first 
chaos transition. All the Smale horseshoes are in the region 0.735 ^ A^ 1.2835. 
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Figure 6: Details of the unstable manifold for A = 0.70. 
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Figure 7: Details of the unstable manifold for A = 0.70. 




Figure 8: Unstable manifold for A = 0.71. The SFP is marked with a black 
circle. 
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Figure 9: Details of the unstable manifold for A = 0.71. 
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Figure 10: Unstable manifold for A = 0.735. The SFP is marked with a black 
circle. 
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Figure 11: Unstable manifolds for A = 0.74, A = 0.75, and A ^ 0.85. The SFP 
is marked with a black circle. For A = 0.85, the four attracting periodic points 
(marked with circles) are near the main fold zones. 
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Figure 12: Stable (black) and unstable (blue) manifolds for A 
A = 0.735, A = 0.74. The SFP is marked with a black circle. 



0.70, A = 0.73, 



These horseshoes are "quadruple" since each branch of each invariant manifold 
intersects the two branches of the other manifold. 

For A = 0.74 the first set of horseshoes has been already created (Fig. 12), 
and for A w 0.748 occurs the sudden expansion of the size of the attractor, 
causing new intersections between the stable and unstable manifolds and new 
horseshoes (Fig. 



13 



A = 0.75). From this point, chaotic and non-chaotic 
zones alternate, because homoclinic tangencies between the stable and unstable 
manifolds are not continuously formed. For example, for A « 0.785 the chaos 
disappears because there are no homoclinic tangencies (Fig. 13 A = 0.80). 
This non-chaotic zone (Fig. 13 A — 0.85, A = 1.05) remains until A sa 1.11 
(Fig. 14 j4 = 1.10, j4 = 1-15), when new homoclinic tangencies are formed. 
This alternation of chaotic and nonchaotic zones is maintained until A ss 1.2835 



(Fig. lA A — 1.21 for a non-chaotic zone and A = 1.25 for a chaotic zone). 



For A « 1.2835 the last set of horseshoes is undone (Fig. 15), in a process 
inverse to the formation of horseshoes for A « 0.735. In fact, for these two 
values of A, the invariant manifolds are qualitatively similars. 

Finally, for 1.2835 ^ A ^ 1.72 there is not any horseshoe. Moreover, there 
are not homoclinic tangencies, and therefore there is no chaos. 



10 




-2 -1 


1 2 
A=0.85 


0.6 


■■"'■'»''-*S.;:-.,^^ 


S;V '; 


0.4 


s 


5 1' J - 


0.2 


■,_ ( 


7 ' _-— "^"""""^ " 





"wrfc^ 


'/-—-^^"'''''''^ 


0.2 


\ff/ 


""■■" "-'V ;•..:,., ,^.„ ■ 


0.4 


p( 


- 


0.6 


r 


- 




- 


2 -1 


1 2 
A=1,05 


0.6 




/ 


04 


\ 


/ ..---^ ' 


0.2 


~\X 


Jr-^-—"'''''''^''^ 





\n 


/^ ' 


-0.2 


v/ 


/ 


-0.4 


^ 




-0.6 




■ 



Figure 13: Stable (black) and unstable (blue) manifolds for A — 0.75, A — 0.80, 
A = 0.85, A = 1.05. The SFP is marked with a black circle. 
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Figure 14: Stable (black) and unstable (blue) manifolds for A = 1.10, ^ = 1.15, 
A = 1.21, A = 1.25. The SFP is marked with a black circle. 
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Figure 15: Stable (black) and unstable (blue) manifolds for A — 1.28, A 
1.2835, A = 1.29. The SFP is marked with a black circle. 



5 Algorithms 

• The Poincare map, /, is computed solving the system (II]) using a Runge- 
Kutta method of order 4-5, with a local error of 10~^^. 

• The SFP is obtained with an absolute error of 10^^^, using a method, 
that we shall describe below, based on the fact that the iteration of a 
neighbourhood of the SFP is also a neighbourhood of the SFP. 

This algorithm is as follows: first consider an initial mesh Mg of a square 
interval /q = [a^o,^^] x [2/0,2/9] containing a unique SFP (x* ,y*). 

Then, given a positive eg, we consider the set 

Ci = {(x,y) e Mo : ||/(x,y) - {x\y*)\\ < eo}. 

Next we repeat this process with another mesh Mi of the square interval 
h = [xl,xl] X [2/1,2/1], where 

x\ := inin{x : {x,y) £ Ci} ; x^ :— max{x : {x,y) G Ci} 
2/1 := mm{y : {x,y) GCi} ; j/i := max{2/ : (x, 2/) € Ci} , 

and for another positive ei < eo- Obviously /i C /o. 

Thus, repeating this process, with a suitable choice of meshes M„ and e„, 
we can obtain a sequence of nested squared intervals Iq Z) Ii Z) I2 ■ ■ ■ such 
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that 

+00 



{x*,y*)=f]ln. 



For example, in our case, we have chosen the initial square interval /g = 
[-1.25,-0.25] X [-0.6,0.4] and an initial mesh Mq of 60 x 60 points for 
every A, with cq — 10^^. In the next iterations we have considered meshes 
of 30 X 30 points, and e„ = lO"""^. We have to note that, in this case, it 
is not necessary to use other methods based on model perturbations, as 
the control methods of OGY or Pyragas (see p, fTUl[TT] V 

• In order to compute the invariant manifolds, first we estimate their tan- 
gential slopes at the SFP. This is done by using the eigenvectors of the 
Jacobian matrix Df of the Poincare map at the SFP. 

• To compute the unstable manifold we iterate a segment with length of 
order 10^'^, centered at the SFP, with the appropriate slope previously 
found. The number of iterations may change with the value of A, but it 
is usually between 12 and 24. 

• For the stable manifold, two complementary methods have been used. 

— First, we have iterated twice a small segment (with the appropiate 
slope), using the inverse of the Poincare map. As was pointed out 
in Section |3.1[ this "inverse method" is only valid to compute the 
branches of the stable manifold until they are out of range. So the 
inital segment must be small enough to keep the iterations within 
range. 

— Next, we find the points of a mesh such that, after n iterations, their 
distance to the SFP is less than a given e. In our case we used n = 8, 
a mesh of 900 x 900 points, and lO'^ < e < 10"^ 

6 Final remarks 

The main results obtained are the following: 

• The existence of Smale horseshoes in the forced BvP oscillator (fTl) has 
been explicitly checked, depending on the amplitude of the external force- 
ment A. In particular, the horseshoes appear for A £ (0.735, 1.2835), 
approximately. 

• The chaotic zones of the bifurcation diagram are related to creation and/or 
destruction of horseshoes, i.e., to the existence of homoclinic tangencies 
between the invariant manifolds. 

• The sudden expansion of the attractor seems to be related to the creation 
of the first set of horseshoes. 
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